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We discuss in detail the coupled-channels approach for the large-angle quasi-elastic scattering in 
massive systems, where many degrees of freedom may be involved in the reaction. We especially 
investigate the effects of single, double and triple phonon excitations on the quasi-elastic scattering 
for 48 Ti, 54 Cr, 56 Fe, 64 Ni and 70 Zn+ 208 Pb systems, for which the experimental cross sections have 
been measured recently. We show that the present coupled-channels calculations well account for 
the overall width of the experimental barrier distribution for these systems. In particular, it is shown 
that the calculations taking into account single quadrupole phonon excitations in 48 Ti and triple 
octupole phonon excitations in 208 Pb reasonably well reproduce the experimental quasi-elastic cross 
section and barrier distribution for the 48 Ti+ 208 Pb reaction. On the other hand, 54 Cr, 56 Fe, 64 Ni and 
70 Zn+ 208 Pb systems seem to require the double quadrupole phonon excitations in the projectiles in 
order to reproduce the experimental data. 



PACS numbers: 24.10.Eq, 25.60.Pj, 25.70.Bc, 27.80.- 
I. INTRODUCTION 

It is now well established that the internal structure of 
colliding nuclei strongly influences heavy-ion collisions at 
energies around the Coulomb barrier. In particular, the 
coupling to the collective excitations (rotation and vibra- 
tional states) in the target and projectile nuclei partic- 
ipating in the reaction significantly enhances the fusion 
cross sections for intermediate mass systems 0, Q • Such 
couplings give rise to a distribution of the Coulomb bar- 
rier 0, [2, y] , which can most easily be visualized for re- 
actions involving a deformed nucleus. In this case, the 
nucleus-nucleus potential depends on the orientation an- 
gle of the deformed nucleus with respect to the beam 
direction. Since the orientation angle distributes isotrop- 
ically at the initial stage of the reaction, so does the po- 
tential barrier. The concept of barrier distribution can 
be extended also to systems with a non-deformed target 
Q, where the distribution originates from the coupling 
between the relative motion and vibrational excitations 
in the colliding nuclei and/or transfer processes. Notice 
that, although this concept is exact only when the exci- 
tation energy is zero, to a good approximation it holds 
also for systems with a non-zero excitation energy 0, H[ • 

In Ref. 6], Rowley et al. have argued that the bar- 
rier distribution can be directly extracted from a mea- 
sured fusion cross section af us (E), by taking the second 
derivative of the product E<j{ us (E) with respect to the 
center-of-mass energy E, that is , D ius = d 2 (E<jf us )/dE 2 . 
This method has stimulated many high precision mea- 
surements of fusion excitation function for medium-heavy 
mass systems @, 0] • The extracted barrier distributions 
have revealed that the concept indeed holds and the bar- 
rier distribution itself provides a powerful tool for in- 
vestigating the effects of channel coupling on heavy-ion 
fusion reactions at sub-barrier energies. It has also been 



shown recently that the concept of barrier distribution 

is still valid even for relatively heavy systems, such as 
ioo Mo+ ioo Mo g 

A similar barrier distribution can also be extracted 
from quasi-elastic scattering (a sum of elastic, inelas- 
tic and transfer processes) at backward angles d, 1C|, 
that is a good counterpart of the fusion reaction 11 1. 
In this case, the barrier distribution is defined as the 
first derivative of the ratio of quasi-elastic to the Ruther- 
ford cross sections da qc i/do-R, with respect to E, i.e., 
D qcl = —d(d<T qc \/d<jfi)/dE. Since the fusion and the 
quasi-elastic scattering is related to each other because 
of the flux conservation, a similar information can be ob- 
tained from those processes and the similarity between 
the two representations for barrier distribution has been 
shown to hold for several intermediate mass systems 

Recently, the quasi-elastic barrier distribution has been 
exploited to investigate the entrance channel dynamics 
for fusion reactions to synthesize super-heavy elements 
HH EH EE EH- It has been shown that the concept 
of barrier distribution remains valid even for such very 
heavy systems once the deep-inelastic cross sections are 
properly taken into account. As is expected, the strong 
channel coupling effects on the barrier distribution have 
been observed. 

In this paper, we carry out a detailed coupled-channels 
analysis for large-angle quasi-elastic scattering data for 
48 Ti, 54 Cr, 56 Fe, 64 Ni and 70 Z 
super-heavy elements Z = 

respectively [H, Ell- We especially study the role of 
multi-phonon excitations of the target and projectile nu- 
clei, which has been shown to play an important role in 
quasi-elastic scattering for the 86 Kr+ 208 Pb system [l4j |. 

The paper is organized as follows. We briefly explain 
the coupled-channels formalism for quasi-elastic scatter- 



.+ 208 Pb systems leading to 
104, 106, 108, 110, and 112, 
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ing in Sec. II. We present the results of our systematic 
analysis in Sec. III. We then summarize the paper in Sec. 
IV. 



II. COUPLED-CHANNELS FORMALISM FOR 
LARGE ANGLE QUASI-ELASTIC SCATTERING 

In this section, we briefly describe the coupled- channels 
formalism for large angle quasi-elastic scattering which 
includes the effects of the vibrational excitations of the 
colliding nuclei. The total Hamiltonian of the system is 
assumed to be 



H 



-V 2 + ^ 0) (r) + 



Z P Z T e 2 



2ju r 
~H cxct + V coup (r, £ P , £ T ), 



(1) 



where r is the coordinate of the relative motion between 
the target and the projectile nuclei, /i is the reduced 
mass and £p and £p represent the coordinate of the vi- 
bration in the target and the projectile nuclei, respec- 
tively. Zp and Zt are the atomic number of the pro- 
jectile and the target, respectively, and is the bare 
nuclear potential, which we assume to have a Woods- 
Saxon shape. It consists of the real and imaginary parts, 
Vff(r) = Vb(r) + iWo(r). H cxct describes the excita- 
tion spectra of the target and projectile nuclei, while 
V coup (r, £p, £t) is the potential for the coupling between 
the relative motion and the vibrational motions of the 
target and projectile nuclei. 

In the iso-centrifugal approximation p| [l?], HH > where 
the the angular momentum of the relative motion in each 
channel is replaced with the total angular momentum J 
(in the literature, this approximation is also referred to 
as the rotating frame approximation or the no-Coriolis 
approximation), the coupled-channels equations derived 
from the Hamiltonian (fTJ) read 



J(J + l)h 2 



+4 0) w + 



Z P Z T e 2 



-E 
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u n(r) + ^ Vnn ' ( r ) U n' ( r ) 



(2) 



where e„ is the eigen-value of the operator -ff C xct for the 
n-th channel. V nn > (r) is the matrix elements for the cou- 
pling potential V coup . 

In the calculations presented below, we use the method 
of the computer code CCFULL [TtJ and replace the vibra- 
tional coordinates £p and £p in the coupling potential 
V conp with the dynamical excitation operators Op and 
Ot- The coupling potential is then represented as 

V coup {r, Op, 6 T ) = V c (r, P , 6 T ) + V N (r, 6 P , T ), (3) 



V c (r, Pl O T ) 



V N {r, Op,O t ) 
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(4) 



1 + cxp 



[r-R a -(RpOp+R T OT)\ 



(5) 



Here, Ap and Ap denote the multipolarity of the vibra- 
tions in the projectile and the target nuclei, respectively. 
We have subtracted vff\r) in Eq. |(SJ) in order to avoid 
the double counting. 

If we truncate the phonon space up to the triple phonon 
states (that is, n=0,l,2, and 3), the matrix elements of 
the excitation operator O in Eqs. and §5$ are given 
by 
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(6) 



where (3 is the deformation parameter, that can be es- 
timated from a measured electric transition probabil- 
ity from the single phonon state (n=l) to the ground 
state (n=0). We have assumed the harmonic oscillator 
model for the vibrations, where e n in Eq. ^ is given by 
e n = n Too. 

The coupled-channels equations, Eq. @, are solved 
with the scattering boundary condition for u n (r), 



u„{r) 



(r — > oo) 



where S£ is the nuclear S matrix. Hj \kr) and 

are the incoming and the outgoing Coulomb 
wave functions, respectively. The channel wave num- 
ber k n is given by \j2ji (E — e n )/h 2 , and ki — k ni = 
^/2^iE/K 2 . The scattering angular distribution for the 
channel n is then given by [l8l | 



(7) 



(8) 



with 



/n( )= £ e K(*)+ CT ,(£- 



2J+ 1 

47T 



Yj O (0) 



-2iir 
yhi k n 



0# - £„,„*) + fc(d)S n , ni (9) 



where <Jj(E) and fc{&) are the the Coulomb phase shift 
and the Coulomb scattering amplitude, respectively. The 
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differential quasi-elastic cross section is then calculated 
to be 



dn 



Eda n 
~7k7 



(10) 



We will apply this formalism in the next sec- 
tion to analyze the quasi-elastic scattering data of 



48 



Ti, 54 Cr, 56 Fe, 64 Ni, and 70 Zn+ 208 Pb systems. 



III. COMPARISON WITH EXPERIMENTAL 
DATA : EFFECTS OF MULTI-PHONON 
EXCITATIONS 

In this section, we present the results of our de- 
tailed coupled-channels analysis for quasi-elastic scatter- 
ing data of 48 Ti, 54 Cr, 56 Fe, 64 Ni, and 70 Zn+ 208 Pb systems 
15 , fl(| . The calculations are performed with a version 
1!| of the coupled-channels code CCFULL [jjfy. Notice 
that the iso-centrifugal approximation employed in this 
code works well for quasi-elastic scattering at backward 
angles 10]. In the code, the regular boundary condition 
is imposed at the origin instead of the incoming wave 
boundary condition. We discuss the stability of the nu- 
merical calculations for the massive systems in Appendix 
A. 

The surface diffuseness of the real part of the nuclear 
potential is taken to be a = 0.63 fm, as suggested by 
recent studies on deep sub-barrier quasi-elastic and Mott 
scattering [2(| Hi], H2| ; while the radius parameter to be 
ro = 1.22 fm for all the systems. Notice that a simi- 
lar value for a has been used also in the analysis of the 
recent experimental data for quasi-elastic scattering in 
the 86 Kr + 208 Pb system Q. The depth parameter, 
Vo, is adjusted in order to reproduce the experimental 
quasi-elastic cross sections for each system. The opti- 
mum values of the depth parameter and the resultant 
Coulomb barrier height are summarized in Table [U As 
usually done, we use a short range imaginary potential 
with Wo = 30 MeV, r w = 1.0 fm and a w — 0.3 fm to sim- 
ulate the compound nucleus formation. The results are 
insensitive to these parameters as long as the imaginary 
part of the potential is well confined inside the Coulomb 
barrier. The excitation energy and the corresponding de- 
formation parameter for the single phonon excitation in 
each nucleus included in the calculations are given in Ta- 
ble HB The latter quantity is taken from Refs. [HI, Hi]. 
The radius of the target and the projectile are taken to 



be R T = 1.2A 



1/3 



and R P = 1.2 A 1 / 3 , 



respectively, in 



order to be consistent with the deformation parameters 
[23l [2^ | . All the calculations shown below are performed 
at the scattering angle of 8 c . m . — 170°. We plot the 
quasi-elastic cross sections and barrier distributions as a 
function of the effective energy defined by [lj| 



2E 



sin(6>/2) 
1 + sin(0/2)' 



(11) 



TABLE I: The depth parameter for the real part of the nu- 
clear potential for the 48 Ti, 54 Cr, 56 Fe, 64 Ni, and 70 Zn+ 208 Pb 
systems. The radius and the diffuseness parameters are taken 
to be ro = 1.22 fm and a=0.63 fm, respectively, for all the sys- 
tems. The resultant barrier height energy Vb is also listed. 



System 



48 Ti+ m pb 
54 Cr+ 208p b 
56 Fe+ 208p b 
64 Ni+ 208p b 
70 Zn+ 208p b 



Vfr(MeV) 



VB(MeV) 



88.90 
91.70 
92.85 
95.10 
108.2 



190.50 
205.50 
222.50 
236.25 
249.30 



TABLE II: The properties of the single phonon states included 
in the present coupled-channels calculations, hui and /3 are the 
excitation energy and the dynamical deformation parameter, 
respectively. 



Nucleus 
2U8 Pb 

48 Ti 

54 Cr 



r 



hw (MeV) 







b Fe 
4 Ni 
3 Zn 



3" 
2+ 
2+ 
2+ 
2+ 
2+ 



2.614 
0.983 
0.834 
0.846 
1.346 
0.884 



0.110 
0.269 
0.250 
0.239 
0.179 
0.228 



"taken from Ref. 
''taken from Ref. 



which takes into account the centrifugal energy. We 
calculate the quasi-elastic barrier distributions from the 
cross sections in a similar way as the one used to ob- 
tain the experimental barrier distributions (l6l | . Namely, 
we use the point difference formula with the energy step 
of AE—0.25 MeV and then smooth the resultant bar- 
rier distribution with the Gaussian function with the full 
width at half maximum (FWHM) of 1.5 MeV. We have 
checked that the shape of the barrier distribution does 
not change significantly even if we use a larger energy 
step for the point difference formula, e.g. A£'=0.5 MeV. 



A. Effect of double phonon excitations 

Let us first discuss the effect of double octupole phonon 
excitations in the 208 Pb target. Such excitations have 
been shown to play a significant role in the sub-barrier 
fusion reaction between 16 and 208 Pb nuclei (2o1. [2o|. 

The dotted line in Fig. Q] shows the result of the 
coupled-channels calculation for the 48 Ti + 208 Pb system 
obtained by taking into account the coupling to the single 
octupole phonon state in the target nucleus, 208 Pb, and 
the single quadrupole phonon state in the projectile nu- 
cleus, 48 Ti. The mutual excitations in the projectile and 
the target nuclei are fully taken into account in this cal- 
culation as well as in all the other calculations presented 
in this paper. Figs. Ufa) andQJb) show the ratio of the 
quasi-elastic to the Rutherford cross sections, da qc i/dan, 
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175 180 185 190 195 200 205 
E eff [MeV] 

FIG. 1: Effects of multi-phonon excitations on the quasi- 
elastic scattering cross section (upper panel) and on the quasi- 
elastic barrier distribution (lower panel) for the 48 Ti+ 208 Pb 
system. The dotted line is the result of the coupled-channels 
calculations including coupling to the one quadrupole phonon 
state in the projectile and the one octupole phonon states in 
the target nucleus, while the solid line is obtained by includ- 
ing the coupling in addition to the two octupole phonon state 
in the target nucleus. The dashed line is the result of double 
quadrupole phonon couplings in the projectile and the double 
octupole phonon couplings in the target nucleus. The exper- 
imental data are taken from Ref. fl6( ] . 

and the quasi-elastic barrier distribution, D qcl , respec- 
tively. Although the overall width of the barrier distribu- 
tion is reproduced reasonably well with this calculation, 
the detailed structure is somewhat inconsistent with the 
experimental data. The situation is similar even when 
we include the double quadrupole phonon state in the 
projectile while keeping the single octupole phonon cou- 
pling in the target nucleus (not shown). We then in- 
vestigate the effect of the double octupole phonon cou- 
plings in the target nucleus. The solid and the dashed 
lines in Fig. Q] show the results with the single and the 
double quadrupole phonon excitations in the projectile, 
respectively, where as the double octupole phonon cou- 
pling in the target is included in both the calculations. 
The former calculation reproduces both the cross sections 



TABLE III: The value of x 2 for the quasi-elastic cross sec- 
tions for the 48 Ti+ 208 Pb system obtained with the coupled- 
channels calculations with various coupling schemes. The cou- 
pling schemes are denoted as [ri2,n3], where is the num- 
ber of quadrupole phonon excitation in the projectile nucleus 
while nz the number of octupole phonon in the target nucleus. 



System 


[1,1] 


[1,2] 


[1,3] 


[2,2] 


48 Ti+ 2l,8 Pb 


19.15 


9.82 


7.12 


37.12 



and the barrier distribution reasonably well, although the 
latter calculation somehow worsens the agreement. This 
clearly suggests that the double octupole phonon exci- 
tations in the target nucleus is important in the quasi- 
elastic 48 Ti+ 208 Pb scattering. We summarize the x 2 
value of our calculations in Tabic IIIII 

Since the coupling to the one quadrupole phonon state 
in the projectile and the two octupole phonon states 
in the target reasonably well reproduce the experimen- 
tal quasi-elastic scattering data for the 48 Ti-|- 208 Pb sys- 
tem, one may expect that the same coupling scheme 
accounts for the experimental data for the other sys- 
tems, 54 Cr, 56 Fe, 64 Ni and 70 Zn+ 208 Pb. The results of the 
coupled-channels calculations with this coupling scheme 
is shown by the dashed line in Fig. [2] Figs. Ufa), (c), (e), 
and (g) are for the quasi-elastic cross sections for the 
54 Cr, 56 Fe, 64 Ni and 70 Zn+ 208 Pb systems, respectively, 
while Figs.[2jb), (d), (f), and (h) are for the quasi-elastic 
barrier distributions. One can clearly see that these cal- 
culations underestimate the experimental cross sections 
at high energies, although the experimental barrier distri- 
butions themselves are reproduced reasonably well. We 
repeat the same calculations by including the coupling 
up to the double quadrupole phonon states in the pro- 
jectile, in addition to the double octupole phonon states 
in the target nucleus. These results are shown by the 
solid in Fig. O The agreement with the experimental 
data is considerably improved, especially for the quasi- 
elastic cross sections. See Table HVl for the \ 2 values. It is 
thus evident that the coupling to the double quadrupole 
phonon states in the projectile is needed in order to ex- 
plain the experimental data for the 54 Cr, 56 Fe, 64 Ni and 
70 Zn+ 208p b reac tions. 

The reason why the double quadrupole phonon cou- 
pling is not necessary for the Ti projectile while it is 
for the heavier projectiles is not clear at this moment. 
This might reflect some ambiguity of the Monte Carlo 
reaction simulation code LINDA [27| which was used to 
subtract the deep-inelastic component from the exper- 
imental yields at backward angles jlq |. Clearly, a fur- 
ther investigation is still necessary concerning the effect 
of deep inelastic scattering on quasi-elastic scattering in 
massive systems [HI, [H], UM, [H| • 
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Cr(1ph)-Pb(2ph) 
Cr(2ph)-Pb(2ph) 
exp. 



54 Cr + 208 Pb 



H h 



H 1 1 1 h 



Fe(1ph)-Pb(2ph) 
Fe(2ph)-Pb(2ph) 
exp. 




56 Fe + 208 Pb 




190 195 200 205 210 215 220 210 215 220 225 230 235 240 
E eff [MeV] E eff [MeV] 



Ni(1ph)-Pb(2ph) 
Ni(2ph)-Pb(2ph) 

Jiii ex P 



H 1 h 



Zn(1ph)-Pb(2ph) 
Zn(2ph)-Pb(2ph) 
exp. 




70 Zn + 208 Pb 




220 225 230 235 240 245 250 255 235 240 245 250 255 260 265 
Eeff t MeV l E eff [MeV] 



FIG. 2: The quasi-elastic scattering cross sections ((a), (c), (e), and (g)) and the quasi-elastic barrier distributions ((b), (d), 
(f), and (h)) for the 54 Cr, 56 Fe, 64 Ni and 70 Zn+ 208 Pb systems obtained with two coupling schemes as indicated in the insets. The 
dashed line is obtained by including the one quadrupole phonon state in the projectile nuclei while the solid line is obtained with 
the double phonon couplings. The double octupole phonon excitations in the target nucleus is included in all the calculations. 
The experimental data are taken from Ref. [la ] . 
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175 180 185 190 195 200 205 
E eff [MeV] 

FIG. 3: Effects of triple phonon excitations on the quasi- 
elastic scattering cross section (upper panel) and on the quasi- 
elastic barrier distribution (lower panel) for the 48 Ti+ 208 Pb 
system. The dashed line is the result of the coupled-channels 
calculations taking into account the coupling to the one 
phonon state in the projectile and the two phonon states in 
the target nuclei. The solid line is obtained by including the 
coupling to the one phonon state in the projectile and the 
three phonon states in the target. The experimental data are 
taken from Ref. [la |. 



TABLE IV: Same as Table UTTl but for the 54 Cr, 56 Fe, 64 Ni and 
™Zn+ 208 Pb systems. 



System 


[1,2] 


[1,3] 


[2,2] 


[2,3] 


b4 Cr+ m Pb 


52.47 


49.80 


20.61 


11.78 


56 Fe+ 20Sp b 


28.46 


28.36 


10.44 


10.28 


64 Ni+ 208p b 


57.45 


61.43 


32.21 


30.64 


70 Zn+ 208p b 


26.52 


24.81 


11.36 


6.87 


B. 


Effect of triple 


phonon 


excitations 





In the previous subsection, we have shown that the 
double octupole phonon excitations in the 208 Pb target 
play an important role in quasi-elastic scattering for the 
systems considered in this paper. However, the calcu- 
lated quasi-elastic barrier distributions have a much more 



prominent peak than the experimental distribution at 
high energies. Since it has been shown in Refs. [HI, [3] 
that the triple octupole phonon excitations of the 208 Pb 
play a significant role in the large-angle quasi-elastic scat- 
tering between 86 Kr and 208 Pb nuclei, it is intriguing to 
investigate such effects in the present systems as well. 

The results of the coupled-channels calculations includ- 
ing the coupling to the triple octupole phonon states in 
208 Pb for the 48 Ti+ 208 Pb reaction is presented in Fig. [3] 
The dashed line is the same as the solid line in Fig.Q] that 
is the result of single phonon in 48 Ti and double phonon 
in 208 Pb. The solid line denotes the results of the triple 
phonon coupling in the target in addition to the single 
phonon in the projectile. By including the triple octupole 
phonons in the target nucleus, the quasi-elastic cross sec- 
tions are improved slightly (see also Table HTT|) . On the 
other hand, one can see that the agreement for the bar- 
rier distribution with the experimental data is much more 
improved by the triple phonon coupling. 

The results for the other systems, 54 Cr, 56 Fe, 64 Ni and 
70 Zn+ 208 Pb reactions, are shown in Fig. |U Figs. [Ha), 
(c), (e) and (g) are for the quasi-elastic cross sections, 
while Figs. [U^b), (d), (f) and (h) for the quasi-elastic 
barrier distributions. Let us first discuss the calculations 
with the single phonon excitation in the projectile. The 
dotted line in the figures is obtained by taking the cou- 
pling to the single phonon state in the projectile and the 
triple octupole phonon excitations in the target. This cal- 
culation underestimates the quasi-elastic cross sections at 
high energies and the obtained barrier distribution is in- 
consistent with the experimental data. Therefore, the 
previous results shown in Fig. [2] is not improved even if 
the triple phonon excitations in the target is taken into 
account as long as only the single phonon excitation is 
considered for the projectile nucleus. The results with 
the double phonon couplings in the projectile together 
with the triple phonon excitations in the target are then 
shown by the solid line in the figure. For comparison, 
we also show by the dashed line the results of the double 
phonon excitations in both the projectile and the target 
nuclei, which is the same as the solid line in Fig. [5] 
One can observe that the inclusion of the triple octupole 
phonon excitations in the 208 Pb somewhat improves the 
agreement between the calculations and the experimen- 
tal data for both the quasi-elastic cross sections and the 
barrier distributions (see also Table [TV)) . 

In Ref. [HI, Ntshangase et al. reduced the coupling 
strength of (3") -> (3 - ) 2 states in 208 Pb by a factor 
of (0.6) and that of (3~) 2 -> (3~) 3 by (0.6) 2 in order 
to explain the experimental barrier distribution for the 
86 Kr-|- 208 Pb reaction. In order to see whether such re- 
duction of the coupling strengths improves the agreement 
between the coupled-channels calculations and the exper- 
imental data for the present systems, we repeat the cal- 
culations by including those effects for the 70 Zn+ 208 Pb 
system. The results are shown in Fig. [51 The solid is ob- 
tained by reducing the coupling strengths as Ntshangase 
et al. did, while the dashed line is the same as the solid 
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FIG. 4: Effects of triple phonon excitations on the quasi-elastic cross sections ((a), (c), (e), and (g)) and on the quasi-elastic 
barrier distributions ((b), (d), (f), and (h)) for the 54 Cr, 56 Fe, 64 Ni and 70 Zn+ 208 Pb systems. The dashed line is the same as the 
solid line in Fig. [2]while the dotted line is the results of the calculations taking the coupling to the triple octupole phonon in the 
target and the one quadrupole phonon state in the projectile nucleus into account. The solid line is obtained by including the 
coupling to the double quadrupole phonon states in the projectile and the triple otcupole phonon states in the target nucleus. 
The experimental data are taken from Ref. [la ]. 
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FIG. 5: Effect of anharmonic octupole phonon excitations 
in 208 Pb on (a) the quasi-elastic cross section and (b) the 
quasi-elastic barrier distribution for the Zn+ Pb reaction. 
The solid line is the coupled-channels calculation obtained by 
reducing the coupling strengths to multi-phonon states, while 
the dashed line denotes the results in the harmonic limit. The 
experimental data are taken from Ref. 



FIG. 6: Comparison of the experimental data with the 
coupled-channels calculations obtained using different val- 
ues of the surface diffuseness of the nuclear potential for 
64 Ni+ 208 Pb reaction for (a) the quasi-elastic cross section and 
(b) the quasi-elastic barrier distribution. The solid and the 
dashed lines is obtained using the surface diffuseness of the 
nuclear potential a = 0.63 fm and a — 1.0 fm, respectively. 
The dotted line is the results obtained by shifting the barrier 
height by around + 2.40 MeV for the calculation using a — 1.0 
fm. Experimental data is taken from Ref. [l(|. 



line in Figs. |4jg) and (h), that is obtained by assuming 
the harmonic limit. In both cases, we take into account 
the coupling to the double quadrupole phonon excitations 
in the projectile nucleus. One can see that the harmonic 
model leads to a better agreement with the experimen- 
tal data both for the cross sections and the barrier dis- 
tribution, as compared to the anharmonic calculation. 
The difference between Ref. [14| and the present calcula- 
tion concerning the role of anharmonicity may originate 
from the fact that Ref. [3| used a smaller value for Rt 

1/3 

(=1.06A T ' fm) and thus a larger value for /3 3 (=0.16). In 
order to clarify the role of anharmonicity of multi-phonon 
excitations in quasi-elastic scattering in massive systems, 
it would be required to take into account also the reori- 
entation terms [HI, [H, . It is beyond the scope of this 
paper, and we will leave it for a future study. 



C. Surface diffuseness of the nuclear potential 

We next discuss the dependence of the quasi-elastic 
scattering on the surface diffuseness parameter of the nu- 
clear potential. The standard value for the diffuseness 
parameter is around 0.63 fm [3(| HH, H2| . Recently, sys- 
tematic studies on quasi-elastic scattering as well as Mott 
scattering at deep sub-barrier energies have revealed that 
the surface region of the nuclear potential is indeed con- 
sistent with the standard value of the surface diffuseness 
parameter [2(| HH [HJ . On the other hand, it has been 
known for some time that the recent high precision data 
of sub-barrier fusion cross sections require a larger value 
of surface diffuseness parameter, ranging between 0.75 
and 1.5 fm [33j. Since the large-angle quasi-elastic scat- 
tering around the Coulomb barrier may probe both the 
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surface region and the inner part of the nuclear potential, 
it is interesting to study the sensitivity of quasi-elastic 
cross sections and barrier distributions to the surface dif- 
fuseness parameter. 

For this purpose, as an example, we repeat the coupled- 
channels calculation for the the 64 Ni+ 208 Pb reaction us- 
ing the nuclear potential with et=1.0 fm. We readjust the 
depth and the radius parameters to be Vq = 160.70 MeV 
and ro = 1.10 fm, respectively, so that the barrier height 
remains the same as the one listed in Table I. We include 
the coupling to the double phonon states in the projec- 
tile and the triple phonon states in the target. Fig. [6] 
compares the results with a=0.63 fm (the solid line) to 
the one with a=1.0 fm (the dashed line). One sees that 
the calculations with a=1.0 fm underestimate the quasi- 
elastic cross section, although the shape of barrier distri- 
bution itself is similar to the one obtained with a = 0.63 
fm. The dotted line is obtained with the same value of 
surface diffuseness parameter a=1.0 fm as the one for the 
dashed line, but by changing the depth parameter Vb so 
that the resultant barrier height is higher by 2.4 MeV. 
This calculation now reproduces the experimental quasi- 
elastic cross sections at energies larger than E c . m , = 230 
MeV reasonably well, but below this energy the cross 
sections are underestimated. Therefore, it seems dif- 
ficult to reproduce the experimental quasi-elastic cross 
sections with the diffuseness parameter of a=1.0 fm at 
energies below and above the Coulomb barrier simulta- 
neously. We have checked that the situation is similar for 
the other systems, 48 Ti, 54 Cr, 56 Fe and 70 Zn+ 208 Pb. This 
result clearly indicates that the standard value of surface 
diffuseness, a = 0.63 fm, is preferred by the experimental 
quasi-elastic data for the systems studied in this paper. 



IV. SUMMARY 

We have performed a detailed coupled- channels anal- 
ysis for large-angle quasi-elastic scattering of the 
48 Ti, 54 Cr, 56 Fe, 64 Ni and 70 Zn+ 208 Pb systems, where 
their experimental barrier distributions have been ex- 
tracted recently. Our coupled-channels calculations with 
multi-phonon excitations in the colliding nuclei repro- 
duce the experimental quasi-elastic cross sections as well 
as the barrier distributions, indicating clearly that the 
coupled-channels approach still works even for massive 
systems [8[ . It was crucial to subtract properly the deep- 
inelastic components from the total backward-angle cross 
sections in order to reach these agreements between the 
calculations and the experimental data. 

In more details, the calculation with the single 
quadrupole phonon excitation in 48 Ti and the triple oc- 
tupole phonon excitations in 208 Pb reproduces reason- 
ably well the experimental data for the 48 Ti+ 208 Pb sys- 
tem. On the other hand, for the 54 Cr, 56 Fe, 64 Ni and 
70 Zn+ 208 Pb systems, we found that the coupling to the 
double quadrupole phonon excisions in the projectile nu- 
cleus in addition to the coupling to the triple octupole 



phonon in the target nucleus seems to be needed to fit the 
experimental data. These results suggest that the triple 
octupole phonon excitations in the 208 Pb nucleus plays 
an important role in describing the experimental data for 
the quasi-elastic cross section and the quasi-elastic bar- 
rier distribution for the present massive systems. This is 
consistent with the previous finding for the 86 Kr + 208 Pb 
system jlil ]. 

Although our calculations well reproduce the gross fea- 
tures of the experimental barrier distributions, higher 
precision data are still required in order to study the 
detailed structure of the barrier distributions, especially 
the role of multi octupole phonon states in 208 Pb. From 
the theoretical side, a further detailed investigation will 
also be necessary, taking into account the anharmonicity 
of the multi-phonon excitations. 
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APPENDIX A: NUMERICAL STABILIZATION 
OF COUPLED-CHANNELS CALCULATION 

In this Appendix, we discuss the problem of numer- 
ical instability of coupled-channels calculations and the 
stabilization methods which we employ in the present 
calculations. 

The coupled-channels equations @ form a set of N 
second order coupled linear differential equations, where 
-/V is the dimension of the coupled-channels equations. 
These equations can be solved by generating N linearly 
independent solutions and taking a linear combination 
of these N solutions so that the asymptotic boundary 
condition, as well as the regular boundary condi- 
tion at the origin, are satisfied. The linearly independent 
solutions can be obtained by taking N different sets of 
initial conditions at r = 0. We denote these solutions 
by ^ m ,(r), where n refers to the channels while rii refers 
to a particular choice of the initial conditions. A simple 
choice for the N initial conditions is to impose 

(f>nm(r) -> cr J+1 S n>m , forr^O, (Al) 

where c is an arbitrary number and J is the total an- 
gular momentum. With these boundary condition, the 
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coupled-channels equations for (j) nm (r) given by 
2/i dr 2 



h 2 d 2 J{J + l)h 2 , N , Z P Z T e 2 



-E + e. 



2[ir 2 



i(r)=0, (A2) 



are solved outwards up to a matching radius i? max - The 
wave functions u n (r) in Eq. ^ are then obtained as 



»(r) 



.4>nm(f), 



(A3) 



where the coefficients C m are determined so that the 
asymptotic boundary condition |7]) is fulfilled. 

In the classical forbidden region, the scattering wave 
functions exponentially damp as the coordinate r de- 
creases. For the smaller energy, the damping is the 
stronger. Therefore, when the excitation energy e„ is 
finite, the absolute value of the wave functions for each 
channel are different by order of magnitude in the clas- 
sical forbidden region, and thus the wave functions tend 
to be dominated by that of the channel which has the 
smallest excitation energy. This easily destroys the lin- 
ear independence of the N numerical solutions 4> nm , and 
causes the numerical instability. This is a serious problem 
especially when the coupling is strong, as in the massive 
systems which we discuss in this paper. 

Several methods have been proposed in order to stabi- 
lize the numerical solution of coupled-channels equations 
[H, [H, [H, H3] • In the present calculations, we stabilize 
the solutions by diagonalizing the wave function matrix 
4>nm a t several points of r in order to recover the linear 
independence. That is, at some radius r s , we compute 
the inverse of the matrix A nm = </>„ m (r s ) , and define the 
new set of wave functions, 



4>nmir) = \6 nh {r) ■ (A 1 ) 



k 



km • 



(A4) 



The new wave functions <f> obey similar coupled-channels 
equations as Eq. (|A2|) . with the boundary conditions 
given by, 

4> nm {r s -h) = y^„,fc(r - h) ■ {A~ l ) km , (A5) 



071771 {j's ) $n 



(A6) 



Here, h is the step for the discretization of the radial co- 
ordinate, r. These coupled-channels equations are solved 
outwards from r s . The solutions <f> can then be con- 
structed as 4> — A-(f>. We impose this stabilization proce- 
dure up to r = 15 fm with an interval of 1 fm. Although 
this method is similar to those in Refs. [36|, l37|, our 
method is much simpler to be implemented. 

This method is sufficient for intermediate heavy sys- 
tems, such as 16 + 144 Sm. For massive systems, how- 
ever, we still encounter a small numerical instability [38| . 
In order to cure this problem, in addition to the stabi- 
lization method (|A4[) . we also adopt two other methods, 
which are used in the computer code FRESCO [39(. That 
is, we introduce two radii, i? m in and i? cu t- i? m in is the ra- 
dius from which the coupled-channels equations (|A2[) are 
solved, i.e., these equations are solved from r = R m in in- 
stead of r = 0, by setting 4> nm (r) = for r < R m i n - Rent 
is a cut-off radius for the coupling matrix, i.e., the off- 
diagonal components of the coupling matrix V nn i (r) are 
set to be zero for r < R cu t- Both the procedures are jus- 
tified when the absorption is strong inside the Coulomb 
barrier, as in heavy-ion systems, and the results are insen- 
sitive to the particular choice of i? m i n and i? cu t as long as 
they are inside the Coulomb barrier. Typically, we take 
Rmin= 6 fm and i? cu t= 10 fm to obtain reasonable results 
for the present systems (notice that the pocket and the 
barrier appear at e.g., 11.3 and 13.2 fm, respectively, for 
the 64 Ni+ 208 Pb system with the nuclear potential given 
in Table I). 
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